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NO BV BOUNDS FOR APPROXIMATE SOLUTIONS TO P-SYSTEM WITH 

GENERAL PRESSURE LAW 

ALBERTO BRESSAN, GENG CHEN, QINGTIAN ZHANG, AND SHENGGUO ZHU 


Abstract. For the p-system with large BV initial data, an assumption introduced in E] by 
Bakhvalov guarantees the global existence of entropy weak solutions with uniformly bounded total 
variation. The present paper provides a partial converse to this result. Whenever Bakhvalov’s 
condition does not hold, we show that there exist front tracking approximate solutions, with uni¬ 
formly positive density, whose total variation becomes arbitrarily large. The construction extends 
the arguments in [3] to a general class of pressure laws. 


1. Introduction 

A satisfactory existence-uniqueness theory is now available for hyperbolic systems of conservation 
laws in one space dimension with small total variation muM- A major remaining open problem 
is whether the total variation remains uniformly bound or can blow up in finite time for large BV 
initial data. Up to now, only few systems of hyperbolic conservation laws are known, where uniform 
BV estimates hold for solutions with large data nails]. On the other hand, examples with finite 
time blowup have been constructed in mm- However, these systems do not come from physical 
models and do not admit a strictly convex entropy. 

In this paper, we focus on the p-system with general pressure law modeling barotropic gas 
dynamics. 

f ut+p{v)a, = 0, , . 

\ vt-u, = 0, ^ ’ 

where T = l/p>0is the specific volume, p > 0 is the density and u is the velocity of the gas. The 
pressure p{v) is a smooth function of v satisfying 

< 0 and pvv > 0. ( 1 . 2 ) 

In [12] . Nishida proved the global BV existence to (II.Ih with large initial data, for 7 -law pressure 
p = v~^ with 7 = 1. On the other hand, in the case 7 = 8 , various front tracking approximate 
solutions were recently constructed in [4], exhibiting finite time blowup of the BV norm. 

For the p-system with general pressure law, in [3|, Bakhvalov extended the global BV existence 
result for isothermal gas dynamics in [T 2 | to any pressure law p[v) satisfying the Bakhvalov’s 
condition 

3pl^ < 2pypyyy forall u> 0 . (1.3) 

In particular, for 7 -law pressure p = v~'^ with 7 > 0, Bakhvalov’s condition holds if and only if 
7 G (0,1]. In [3], more general 2 x 2 systems of conservation laws are also considered. 

We observe that Bakhvalov’s condition determines whether the strength of a shock increases or 
decreases by crossing a shock of the opposite family, as shown in Figure [H The shock strength is 
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here measured by the change of h{v) across the shock, where 

h{v) = / ^/-pv dv 

J V 

is the density part in the Riemann invariants 

s = u + h{v), r = u — h{v). 


(1.4) 



Figure 1. The upper figure shows two interacting shocks in the x-t plane. For this 
interaction, the lower figures show three different cases, in the (u, /i)-plane Left: the strength 
is amplified after the crossing, and Bakhvalov’s condition m is not satisfied. Middle: the 
strength is same after crossing. Right: the strength is reduced after crossing. In the middle 
and the right pictures, Bakhvalov’s condition (O is satisfied. For 7-law pressure with 
7 > 0, the left figure corresponds to 7 > 1, the middle to 7 = 1, and the right to 0 < 7 < 1. 

In the present paper we extend the blowup examples in [4] to the case where the pressure violates 
Together with [3], this indicates that Bakhvalov’s condition (II.dp is necessary and sufficient 
for the BV stability of the front tracking scheme. More precisely, the following result will be proved. 


Theorem. Assume that the pressure p{-) satisfies for every v > 0 but violates (|1.3I) for 
some V > 0. Then there exists a front tracking approximate solution where the density remains 
uniformly positive while the total strength of waves approaches infinity as t ^ 00 . At each wave- 
front interaction the strengths of outgoing waves are the same as in the exact solution. The only 
errors introduced by the front tracking approximation are in the speeds of the wave fronts. 

By suitably modifying the construction given in the last section of [4j , we expect that one could 
also construct an example of front tracking approximation where the BV norm blows up in with 
finite time. The main ideas leading to the blow-up example can be explained as follows. 

If Bakhvalov’s condition (jl.3ll fails for u in a neighborhood of vq, one can construct two small 
approaching shocks such that 
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(i) their left and right states remain in the region where Bakhvalov’s condition fails, and hence 

(ii) calling cri,cJ2 their sizes before the interaction and cr^,cr2 their sizes after the interaction, 
one has 


(Tl = 172 < a'l = (72 . (1.5) 

Next, assume that these small fronts bounce back and forth between two very large shocks (Fig. [ 2 l 
left). After a hrst reflection at the points Ai, A2, two rarefactions are created. When these rarefac¬ 
tion impinge again on the large shocks at Bi,B2, they generate two new shocks. Every time a front 
is reflected by a large shock, the outgoing wave is strictly smaller than the incoming one. However, 
is the shocks 81,82 are very large, the strengths of incoming and reflected fronts are almost the 
same. Thanks to m, by a suitable choice of the shock strengths, we can achieve 


<Tl 


(^2 = 




= 4 '. 


( 1 . 6 ) 


Hence the interaction pattern can be iterated in time. 

If we further increase the strengths of the shocks S'!, 5*2, in (II.6p we would have 

ki| = k2| < kii = |(72|. (1.7) 

To achieve again a periodic pattern, one needs to cancel part of the rarefaction emerging at A2. 
As shown in Fig. O right, this can be done by merging it with a shock of the same family, at the 
interaction point Ai. In the end, this yields an asymmetric, periodic interaction pattern where 

// // 

(7i — (7l , (72 — (72 . 




Figure 2. Solid lines denote shocks while dotted lines denote rarefaction fronts. Left: 
a symmetric periodic interaction pattern, where two small fronts bounce back and forth 
between two large shocks. In the region between the two large shocks, the solution takes 
values in the region where Bakhvalov’s condition fails. Right: an asymmetric interaction 
pattern. Here part of the rarefaction originating from A2 is canceled at Ai by merging with 
a shock of the same family. 
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Next, on top of these periodic patterns we add an infinitesimally small wave front (a compression 
or a rarefaction), as in Fig. [3l If this additional front is initially located at P and has strength e, 
after a complete set of interactions we show that 

(i) For the symmetric interaction pattern the strength of the small front at Q is e' = e + o(e). 

(ii) For the asymmetric interaction pattern the strength of the small front at Q is e' = Ke + o{e), 
for some k > 1. 

As this cycle of interactions is repeated over and over, the infinitesimal front is enlarged by an 
arbitrarily large factor. 

Finally, as in [1], we replace this infinitesimally small front with a train of countably many pairs 
of rarefaction-compression fronts having sizes ±2“^e, with k = 1, 2,.... This yields a front-tracking 
approximate solution satisfying the properties stated in the Theorem (see Fig. [7|). 




Figure 3. If a small front is added on top of the interaction patterns in Fig. [21 after a 
complete set of interactions the strength of this front is (i) almost the same, in case of the 
symmetric pattern on the left, and (ii) strictly larger, in case of the asymmetric pattern on 
the right. 


Remark. We emphasize that our result does NOT imply that the total variation of entropy 
weak solutions to the p-system can become arbitrarily large. Rather, it shows that front tracking 
approximations can be unstable in the BV norm, whenever Bakhvalov’s condition is violated. 
The present construction also shows that for large initial data, uniform a priori bounds on the 
total variation cannot be proved simply by estimating the wave strengths at each interaction. As 
remarked in [4], to establish such BV bounds (if they do indeed hold) it will be essential to use also 
the decay of rarefaction waves, due to genuine nonlinearity. 

The paper is organized as follows. In Sections 2 and 3 we study the wave curves and calculate 
wave interactions. In Section 4 we first construct a front tracking approximate solution with a 
periodic interaction pattern. Then, by suitably perturbing this periodic pattern, we give examples 
of front tracking approximate solutions where the BV-norm blows up as t —)• oo. 
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2. Wave curves 


In this section, we introduce basic notation and review the rarefaction, compression, and shock 
curves for dni). We omit some standard calculations for wave curves and refer the reader to 
Chapter 17 in [T3] for details. 

In Lagrangian coordinates, the wave speed for is 

c = ^/-p'{v). 

Integrating the eigenvectors of (ll.ip . one obtains the Riemann invariants s and r: 


s = u + h r = u — h , 


with 


h = h{v) = / y/—Pv dv. 

Jv 

In the case of a smooth solution, s and r satisfy 

St + CSa; = 0 , Ti - era; = 0 . 


( 2 . 1 ) 


( 2 . 2 ) 


This yields the curves for the Rarefaction and Compression simple waves. 
For a shock wave, the Rankine-Hugoniot jump conditions take the form 


a[u] = [p(u)] , (2.3) 

a[v] = -[u]. (2.4) 


where [u] = Ur — ui, etc..., and the subscripts I and r denote the left and right states on the shock 
wave, respectively. Together with the Lax entropy condition, this uniquely determines the shock 
curves. 

The following table summarizes the equations for rarefaction, compression and shock curves. We 
refer the reader to Chapter 17 in [H] for detailed calculations. We use {u, v) and {u, v) to denote 

the left and right states across the wave, respectively. Moreover, we use R, R, C, C S, and S 
to denote the forward or backward (or second or hrst) rarefaction, compression and shock waves 
respectively. 


R : 

u — u = h{v) 

- h{v), 

V > 

R : 

u — u = h{v) 

- h{v), 

V < 

C : 

u — u = h{v) 

- Kv)-, 

V < 

C : 

u — u = h{v) 

- h{v), 

V > 

S : 

u — u = —sj {v — v){p{v) — p{v)), 

V < 

S : 

u — u = —-\/{v — v){p{v) — p{v)), 

V > 


We recall that the combined shock-rarefaction curves have regularity mm- 


3. Wave interactions 

In this section, we calculate the head-on interactions between two shocks and between a shock 
and a rarefaction, respectively. 
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3.1. Preliminaries. Consider the function 


a = a 


pv 

{v,v) = h{v)-h{v) = / y/-pvdv. 

J V 


(3.1) 


Since a„(u,u) > 0, one can recover u as a function of a and v, say, v = v{a,v). We also introduce 
the function 


F{a,v) = 


\J(v{a, - v) {p{v) -p{v{a, u))) 


(3.2) 


We now compute the Taylor expansion of p{v), for v near v. In turn, this can be used to calculate 
the Taylor expansion of F{a,v). 

p{v) = p{v)+p'{v){v-v) + ]-p"{v){v-v)'^ + l;p"'{v){v-v)^ + ^p^‘^\v){v-v)'^ + o{v-v)‘^ . (3.3) 

2 b 24 

Using (13.ip and considering v = v{a,v), we compute 

v — v = v(a,v) — v(0,v) (3.4) 

= (-P'(v))~^a + ^(-p'(v))~^p"(v)a^ + i (p"(v)f - ^p'(v)p"'(v) (-p'(v))~h^ + O(a^). 

Using (|3.3|) and (|3.4I) . we obtain 

F(a, v) = \v — u| 


p{v)-p{v) 


v — v 


where 


= sign(u — u) • a {l + Ji{v)a^ + J 2 (u)a^} + o(a^), 

J2 = 


(3.5) 

(3.6) 

(3.7) 


At 



Figure 4. Head-on interactions 

3.2. Head-on wave interactions. In this section, we consider interactions between two opposite 
waves, as shown in Figured! where the incoming waves can be rarefaction, shock, or compression 
waves. The wave does not change its type after crossing a wave of the opposite family. 

We use subscripts 1, 2, 1' and 2' to denote the incoming and outgoing waves of the first and 
second family, respectively. And we denote the {u, v) states between these waves according to 
Figured! For any wave-front, we denote by use a = /ijeft — bright the difference between the values 
of h at the left and right states of the front. For example, referring to Figured! one has 

a2 = h — ho, ai = ho — hj., ai' = h — , 02 ' = ~ hr ■ (3.8) 
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Shock-Shock interaction. In this part, we consider the weak shock-shock interaction. For sim¬ 
plicity, we only consider the case when two shocks have same strength, i.e. v = Vr, hence h = hr- 
Using (12.5p and ()3.8p with h = hj., one obtains 


0-2 
02 ' 

Then by (12.5p and (13.2p , we have 
which yields 

02 ' = 

where J 2 was dehned at (j3.7l) . 


h — ho = hr — ho = —ai > 0 , 
hm hr — hm h — (X\' 0 

F{a 2 ,v) = F{a[,v) 

-ai' = 02(^1+ 2 J 2 {v)afj + 0(02), 


(3.9) 


(3.10) 


3.3. Rarefaction-Shock or Compression-Shock interaction. We now consider the interaction 
between a backward rarefaction and a forward shock. By (|2.5I) . we know ai, 02 , oi', and 02 ' are all 
negative. Traversing the waves before and after interaction yields 


—aii+F(ai',v] = —oi+F 


ai',VQ{a 2 ,v)y 


By the equation (13.51) . we thus have 

ov = ai + J 2 “ 2 a? + 0(020^). 


Hence, by 


01+02 = Oi' -|- 02 ' 


we obtain 

I 02 I = |o2|(l +-^2|oi|^) + o(|a2af|) (3T1) 

By an entirely similar calculation, we have same estimate for the interaction between a backward 
compression and a forward shock. By symmetry, a similar estimate holds for the interaction between 
a forward compression and a backward shock. 


4. Front tracking approximations with unbounded BV norm 

In this section, we construct a front tracking approximate solution whose BV norm tends to 
infinity as t —)• 00 . We assume that the Bakhvalov condition ()1.3p fails at some u > 0. Hence, by 
continuity and by (jl.2p there exists some interval {vl,vu), in which 

J2 = -^P” {^pIv - ‘^PvPvvv) > 0 foT all V G {vL,vu)- (4.1) 

4.1. Front tracking approximations with a periodic interaction pattern. Following [3], 
we hrst construct a symmetric interaction pattern containing four wave fronts, as shown in Fig. [5j 
This pattern is symmetric, because two boundary shocks Si and S 2 (and also the inner shocks AiC 
and A 2 C) are chosen to have the same strength measured by the difference in h between two sides 
of each shock. We choose the strengths of the two large shocks ^i, S 2 and of the two intermediate 
waves in such a way that, after a whole round of interactions, these strengths are the same as at 
the initial time. Working in the (tt, h) plane, this is achieved as follows. 

(i) Choose states Ai, A 2 , Bi, B 2 , C and D such that v G {vl,vu) at these states. Hence (|4.1I1 
is satished inside and on a neighborhood of the diamond with vertices Hi, C, H 2 , D. 

(ii) Construct two shocks: the 1-shock AiC and the 2-shock A 2 C, approaching each other. 
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Figure 5 . A periodic interaction pattern. The left picture is on h-u plane. The right 
picture is on t-x plane 


(iii) Determine the two outgoing shocks DAi and DA 2 , resulting from the crossing of the above 
two shocks. 

(iv) Construct a rectangle having two opposite vertices at C and D. Call Bi, B 2 the remaining 
two vertices. 

(v) Finally, the state Ui is chosen so that the two points Bi and A 2 are on the same 1-shock 
curve with left state C/p Symmetrically, Ur is chosen so that the two points B 2 and Ai are 
on the same 2-shock curve with right state Ur- 

We observe that, by (j.S.lOjl and (j4.ip . the /i-component of the states Bi and B 2 is larger than the 
/i-component of Ai and A 2 . 

The existence of states Ui, Ur satisfying (v) is now proved in the following lemma, illustrated in 
Fig. [6l 

Lemma 4.1. In the {u, h)-plane, consider two points Bi = {ui,hi) and A 2 = (u 2 ,/i 2 ). Assume 
that 

(i) ui < U 2 , and hi > h 2 . 

(ii) Calling A = ( 112 )^ 2 ) the point on the 1-shock curve with right state Bi with the same 
u-component as A 2 , one has < h 2 - 

Then there exists a unique Ui = {ui,hi), with 0 < hi < /i 2 , such that both Bi and A 2 lie on the 
1 -shock curve with left state state Ui. 

Remark 4.2. Condition (ii) clearly holds when the interaction diamond A 1 -C-A 2 -D is small 
enough, i.e. the interactions inside the diamond are all weak. 

Proof. We shall use (13.21) with {ui,vi) while {u,v) = (tti,ui) or ( 1 x 2 , 112 ). To prove the lemma we 
need to hnd {ui,pi) such that 

ui-ui = {p{vi) - p{vi)) {vi - ui), U 1 -U 2 = \J{p{v 2 ) - p{vi)) {vi - V 2 ). (4.2) 

This will be achieved if we can find vi such that 

U 2 -U 1 = G{vi ), 

where G is the function defined as 

G{v) = (p(ui) - p{v)) [v-vi)- {p{v 2 ) - p{v)) {v - V2) ■ 


( 4 . 3 ) 
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The assumption (ii) implies 

G{v2) = (p(ni) - P{V2)) {V2 - ui) < U 2 -U 1 . 
Moreover, a direct computation shows 


lim G{vi) = + 00 . 

Vl^OO 


Finally, for any < ^2 < n, we have 


9 _ -p'{v){v - vi) + p{vi) - p{v) -p'{v)(v - V 2 ) + p{v 2 ) - piv) 

—G{v) - ^ > °' 


2v/ {p{vi) - P{v)) (-y - III) 2 J {p{v2) - p{v)) {v - V2) 


Indeed, since p''{v) > 0, one has 




. _ , - 3 <0 

2 ^{p{a)-p{v)){v-a) ^ {p{a) - p{v)) {v - af 

for any a < v. Since vi > V 2 , there exists a unique value of vi such that G{vi) = U 2 — ui. 


□ 



Figure 6 . By moving the point U along the 1 -shock curve with right state Si, we even¬ 
tually reach a left state Ui such that the 1-shock curve through Ui contains A2 as well. 



Figure 7 . A periodic pattern that amplifies a train of small wave fronts. 























10 


ALBERTO BRESSAN, GENG CHEN, QINGTIAN ZHANG, AND SHENGGUO ZHU 


4.2. An example with unbounded BV-norm. Next, as shown in Fig. [71 on top of the periodic 
pattern constructed in Fig.[5l we add a train of countably many pairs of rarefaction and compression 
waves. The /c-th pair of waves have sizes ±2“^e. Notice that if a front of arbitrary size a crosses a 
rarefaction and then a compression wave of exactly opposite sizes, after the two crossings the size 
of the front is still u, exactly as before (Fig. [71 center). As a result, the interaction pattern of four 
large fronts retains its periodicity. 

Note that in Fig. [71 we perturb the symmetric periodic pattern in Fig. [5l to an asymmetric 
periodic pattern by splitting some reflecting rarefaction wave into two pieces. The detail of this 
perturbation will be discussed later. We recall that the strength of a wave is always defined as 

|Q(^left)'bright) I I bright ^leftl ) 

where the subscripts denote the left and right states across the wave-front, respectively. 

We always assume that each front in the train of small waves has strength < e. Indeed, we can 
always perform a partial cancellation of the compression-rarefaction pair so that both fronts have 
strength < e. We choose e > 0 small enough so that all states between two boundary shocks satisfy 
V G ivL,vu)- 

We consider the amplification of total wave strength of these alternating waves. To fix the ideas, 
consider a 1-rarefaction or compression of strength sa > 0, located at A. Within a time period, 
this front will 

i. Cross the intermediate 2-shock. 

ii. Interact with the large 1-shock at Pi producing a 2-compression. 

iii. Cross the intermediate 1-shock. 

iv. Cross the intermediate 1-rarefaction. 

V. Interact with the large 2-shock at P 2 producing a 2-rarefaction. 

vi. Cross the intermediate 2-rarefaction. 

Indeed, when a small wave of strength e~ crosses a shock of the opposite family of strength s, 
by (13.lip the strength of the outgoing front is 

= ( 1 -I-J 2 S^-I-o(s^)) e“. (4.4) 

When the front crosses a rarefaction of the opposite family, its strength does not change. 

Finally, when the small wave impinges on a large shock at Pi or at P 2 , we need to estimate the 
relative size of the reflected wave front. 

Calling £~, the strengths of the front before and after interaction, to leading order we have 

= (1 —2tan0)e“ (4.5) 

where 0 is the angle between line segments B 1 A 2 and A 1 A 2 in Figure [5l and s is the strength of 
inner shocks AiC or A 2 C. 

When the additional front reaches P, we want its size to be increased by a factor k > 1. To 
achieve this goal, we need to perturb the symmetric periodic pattern into an asymmetric periodic 
pattern as shown in Figure [8| 

As in the figure, for simplicity we assume he = l,ttc = 0, = r, ui, = s. Using the Rankine- 

Hugoniot condition, we can calculate and h^. 

Indeed, Ub - Uc = s = ■\/-[p(r’fe) - p{vc)]ivb - Vc). 

s = \l-\p'{vb - Vc) + ^p"{vb - VcY -h - Vef + ^P^'^^Vb - Uc)^ 0{{vb - Vc)^)]{vb - Vc) 

(4.6) 
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Figure 8 . Amplification of infinitesimal waves. 


By expressing Vb — Vc in powers of s, one obtains 


Vb 


-Vc = l + \/^s + -^s^ + 




//2 


+ 


1 


4^2 l32(-p')V2 12 (-p')5/2 

Considering (jl.4p . we have 




^//3 p^^'^ 


8 p'^ 8 p'^ 48p'^ 




hb — 1 — — s + 
In a similar way, we obtain 
hd-I = r - 


1 


+ 


p 


//3 


+ 


p//^(3) 


96 {—p'Y \ 64(—p')^/2 95(_p/)7/2 




(4.7) 


1 p "2 


:r^ + 


p 


//3 


+ 




+ C>(r®) 


96 (-pO^ \64(-p')^/^ 96(-p')^/^ 

Since g is the intersection point of two rarefactions, we can calculate the coordinate of g as 


(4.8) 


1 p"2 


Ug = {s + r)-—j^{s +r)+0(r,s). 


ho = l + r — s + 

® 192 {-p') 


,112 Y 


1 P‘ /A 3 


-O + x 


p 


//3 


+ 


p"p(^) 


(^4 ^g4)^ 0(^5^ ^5)^ 


5 


2 y64(—p')®/2 96(—p')'^/2 

Hence the slope of eg is 

r- s + - r^) + i (g4(i(p,)9/2 + 96[_p/)7/2) (t""^ + s^) + ^>(5^, t-®) 


tan 6 = 


S + r - 29 P"\, (g 3 _L ^ 3 ') 

* ^ ' 192 (-p')a I* w I ) 


(4.9) 


(4.10) 
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For the left boundary shock, we can repeat above process. By assuming Ua — = s, Ua — = r, 

we can obtain the slope of ae, which is similar to M.lOh . 


ton 9 = 


r — s + 


192 


(^3 _^3) U 


r,(3) 


P P 

2 ' 


)(^4 _|_ g4^ _|_ 0(^55^ 




(4.11) 


Here r and s are independent, so we can take different relations between r(r) and s{s) for the right 
and left shocks. From the figure, we expect r < s. 

For the right boundary shock, we take r = s — ( g 4 (^p/) 9/2 + 9 g^_p/) 7/2 To leading order, the 
slope of the shock curve eg is 

tan0 = o(s^). 

We take f = s — 2J2S^, so the slope of ae is 

tan^ = {-J2 + + o{s^). 

Since r — s = f — s, the relation between s and s is 

After one complete set of interaction, the strength of the small wave located at B is 

£B = (1 +J2S^+o(s^))^ ^1 — 2J2S^ + {~p') (l+o(s^))ey4 = (1 + + o(s^)) £a ) 


where 

/ 1 \ 3/4 

X = ((-p')-’/^ - 2)J, + 2J. > 0. 

The small wave has been amplified by a factor 1 + + o(s^). 

By construction, after each period each pair of small compression-rarefaction wavefronts is en¬ 
larged by a factor > A > 1. When a pair grows to size > e, we can perform a partial cancellation 
so that its size remains G [e/2, e]. After this manipulation, we can restrict the specific volume v 
to be in the interval {vl,vu). So the condition (|4.ip always holds in the construction. 

Since the total number of small wave-fronts is inhnite, after several periods a larger and larger 
number of pairs (compression -|- rarefaction) reaches size > e/2. Hence, as f —>■ oo, the total 
variation of this approximate solution grows without bounds. 


Appendix 

Some detailed calculations about the slope of shock curves are given below. 

As in the figure, for simplicity we assume he = l,Uc = 0, Ud = r, Uf, = s. By R-H condition, 
Ub — Uc = s = y/—\p{vb) — pivc)]{vb — Vc)- By doing Taylor expansion, we have 

\j-[p'{Vb - Vc) + ^p"{vb - VcY + - Vef + ^P^'^^Vb - Uc)^ + 0{{vb - Vc)^)]{vb - Vc) 

(4.13) 


s = 
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We want to express Vh — Vc in powers of s. Assume Vh — Vc = 1 + As + Bs^ + Cs^ + Ds'^ + 0{s^), 
then compare the coefficients of s^, s^, s® to determine constants A, B, C, D. 

s2 = -\p'{As + Bs"^ + Cs^ + Ds^ + 0{s^)) + y'{As + Bs^ + Cs^ + Ds^ + 0{s^))‘^ 


+ Ip^^\As + Bs^ + Cs^ + Ds* + 0{s^)f + ^p^^\As + Bs^ + Cs^ + Ds^ + 0{s^)y + 0{s^)] 


•(As + Bs^ + Cs^ + Ds^ + 0(s5)) 


Coefficient for s^: 


Coefficient for s^: 


Coefficient for s^: 

0 


1 = -^V, 

A = ±(-p')-^/^. 


0 = -A(Bp' + 2^V) - ^Ap', 


B = -^A^ = 


4p' 


4p‘ 


n' 


= -A{Cp' + p"AB + ^ ^A^p”) - CAp' 


(4.14) 


= -2ACp' - p"A^B - _ lA^Bp". 


, p"2 p"2 p"2 5p"2 p(3) 

^ 4p'3 gp/2 16p'3 gp/3 16p'3 gp/2 

5 y/2 1 ^(3) 

32(-p/)7/2 12 (-p')5/2- 

Coefficient for s^: 

0 = -A{Dp'+ yB‘^+p''AC+ \p^^hA^B + ^pA)A^) 
-B{Cp' + p” AB + \p^^^ A^) 

-C{Bp' + \A^p’') 

-DAp' 


2ADp' = -2BCp' - ^p"AB^ - ^p”A‘^C - ^p^^^A^B - 


D =- 


BC 3p 


- ~B^ - -^AC -- 


3 p" 


1 p' 


(3) 


A A p' 


4 p' 


3 p' 


-A^B- 


1 p' 


(4) 


48 p' 


-A^ 


p>l3 p"p(3) p(4) 

8p'5 gp/4 48p'3 
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So by definition of and values of A, B,C, D, we have 

_ /-Kc + A 

hh-l = 


pVc _ rVc-\-^^ 

/ \/^dV = - \/^dt 

J Vh J Vr- 


/ Dc + A -I -I 1 -| 

- •=) + 2 [- 4 (-p')"’''V“ - 2 (-P')^‘''P<”l(» - Vcf 

+^i-l(-pr^'v^ - ^(-p'r^'Vp<’> - j(-p')-’'“pV’> - ^(-pO-'^v-’k* - p.)’* 

= -vV^A - l(-pr''^(-p")lA‘ - - i(-p')-‘''V«lA» 

= s' 

+ (-(-p')''"c + i(-p')-‘'VpiB - j[-j(-p')-’'V" - 1(-pT''V‘'1(-pT‘A s' 

+ ('-(-p')‘"d - i(-p')-‘'"(-p")i(S' + 2AC + + i(-p')-‘'V'>|3sl' 

-2^[-j(-p')-''V" - ^(-p')-'''?''?''’ - 5(-pr‘"p'*’isi‘) s' 


B 


1 P"" 3 / P‘ 

- + +( 


,//3 


.n7/2^^ 


96 (-p')^ A4 (-p')^A ^ 96 (- p ') 7/2 

( 4 . 15 ) 

Since the rarefaction curves dg and bg are perpendicular in the Figure El we can solve the 
following system for {ug,hg). 




»„(3) 


--- = (l + --M(^-^ + (64(^ + 96f^ 

1 p "2 


(4.16) 


u, = {s + r)- + r^) + 0{r\ s% 

1 p 


ha = l + r- s + 


JI 2 1 


p 


//3 


+ 


pip) 


192 (-pO^ 2 \^64(-p')®/^ 96(-p')^/^ 

So the slope of the shock curve rg at g is 

p»2 ^ ^ 

T92 PVF'''®'" “ 2V64(4p')9/2 -E 96(p/)7/2 


(r^ + s^) + C)(r^s^) 


(4.17) 


^ ^ r-s + ^,^is^-r^) + 
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